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Monte Carlo simulations, in which the Schrodinger equation is solved at each Monte Carlo sweep, 
are employed to assess the influence of magnetization fluctuations, short-range antiferromagnetic in- 
teractions, disorder, magnetic polaron formation, and spin-Peierls instability on the carrier-mediated 
Ising ferromagnetism in two-dimensional electronic systems. The determined critical temperature 
and hysteresis are affected in a nontrivial way by the antiferromagnetic interactions. The findings 
explain striking experimental results for modulation-doped p-Cdi-zMnzTe quantum wells. 

PACS numbers: 75.50.Pp, 05.10.Ln, 75.30.Et, 78.55.Et 



Over the recent couples of years, Mn-doped III-V and 
II- VI diluted magnetic semiconductors have become an 
important playground for developing our understand- 
ing of carrier- mediated magnetism in solids Q. To a 
large extend, this stems from the fact that in these sys- 
tems the relevant interactions can be tuned by chang- 
ing the carrier and magnetic ion densities as well as by 
imposing strain, confinement, electric field, or illumina- 
tion. Particularly intriguing are properties of magnetic 
quantum wells (QWs) [2,03,0], as a critical dimension- 
ality is two for a number of pertinent phenomena in- 
cluding the stability of ferromagnetism against magne- 
tization fluctuations U Anderson localization [f|, spin- 
wave Peierls-like instability h| , and formation of self- 
trapped magnetic polarons Nevertheless, in the case 
of modulation-doped p-type (Cd,Mn)Te QWs, the tem- 
perature T c at which spontaneous spin splitting of elec- 
tronic levels appears as well as its temperature depen- 
dence 0, IToL fill fl^ l follow predictions of a simple mean- 
field Zener-like model of ferromagnetism Q • At the same 
time, however, wide hysteresis loops that are expected 
within this model 0] have not been observed. Instead, ac- 
cording to polarization-resolved photoluminescence (PL) 
measurements, the global spin polarization of the carrier 
liquid increases slowly with the external magnetic field 
along the easy axis, reaching saturation at a field by a 
factor of twenty greater than what could be accounted 
for by demagnetization effects flfl fTlj| . 

In order to determine the importance of various phe- 
nomena that control magnetism in such a reduced dimen- 
sionality disordered system, we have undertaken mas- 
sive Monte Carlo (MC) simulations involving solution 
of the Schrodinger equation at each MC sweep. Our 
system consists of a p-doped Cdo.96Mno.04Te QW in 
which a magnetic order has been observed below 3 K 
H OH El 03 The QW of width L w is sandwiched 
between non-magnetic (Cd,Mg,Zn)Te barriers containing 
acceptors that supply holes to the QW. In this modula- 
tion doping arrangement the spatial separation of the 



holes and the acceptors reduces significantly the impor- 
tance of ionized impurity scattering. Moreover, accord- 
ing to effective mass theory and optical studies of this 
confined system @, 0, 0, 0] the in-plane and per- 
pendicular hole motions are decoupled, so that the hole 
wave function can be factorized as ^(i?) = ^p^{r)ip{z), 
where f{z) — y/2/Lw sxd.(ttz/Lw). At the same time, 
the combined effect of confinement and biaxial epitax- 
ial strain breaks rotational invariance, separating by 15 
meV the ground state heavy hole subband L z = ±1 
from the light hole subband L, = [TlL This, to- 
gether with a strong spin-orbit interaction, XsL, where 
A « 0.7 eV, fixes the hole spin along the z axis. Ac- 
cordingly, the coupling between the hole spin and Mn 
spin assumes an Ising- like form 0, 0, Q, 0, 
H p d = —I p d(R— Ri)s z 'Sf , where I p d is the p-d exchange 
function. This strong magnetic anisotropy renders phe- 
nomena associated with quantum fluctuations of the hole 
spins unimportant. 

Within the standard Ruderman-Kittel-Kasuya-Yosida 
(RKKY) approach, the carrier-mediated exchange inter- 
action is evaluated within the linear response theory, 
which is valid as long as the exchange induced spin- 
splitting is small compared to the Fermi energy, and the 
carrier wave functions are not perturbed by the Mn spins. 
Here, we go beyond the RKKY model, as we obtain the 
carrier wave functions and energies for a particular Mn 
spin configuration by exact diagonalization. If the p- 
d exchange energy \/3N \ = \(u k \I P d\uk}\, where u k is 
the Bloch amplitude normalized in the unit cell volume 
1/N — dg/4, is greater than the width of the carrier 
band W w YlhT /m*a 2 , the 3D potential well associ- 
ated with the single Mn atom can bind the carrier [l^j . 
Since in the case of p-Cdi-^Mn^Te, (3N Q = -0.88 eV, 
to* = 0.25to o , a = 0.647 nm, and thus W = 8.8 eV, 
we are in the weak coupling regime, where a single Mn 
spin affects weakly the orbital part of the hole wave 
function. Hence the p-d interaction can be cast in the 
Fermi contact form. Thus, the Hamiltonian describing 
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FIG. 1: Temperature dependence of magnetization and sus- 
ceptibility of holes and Mn ions, for selected values of hole con- 
centrations in a Cdo.96Mno.04Te QW. Open circles : Jy = 0, 
closed circles : Jy / 0. 



the in-plane motion of holes for a given distribution of 
Mn spins is H KE = p 2 /2m* - E, (3S(r - n)\ V (zi)\ 2 s z Sf , 
where Si is a classical spin vector with S — 5/2, and 
s z is the hole spin density operator with s — 1/2. We 
allow for antiferromagnetic (AF) interactions between 
Mn spins, H AF = V, , k !l J,,S,S,. with Jy = -6.3, 
— 1.9, and —0.4 K for the first, second, and third neigh- 
bors, respectively 0]. The total Hamiltonian reads 
H = Hke + Haf- 

In our simulations we consider a finite cell Lx Lx Lw, 
cut from the (Cd,Mn)Te FCC lattice. We note that since 
the hole densities in the system in question are up to 
10 4 lower than that of the Mn ions we have to employ 
much larger cells than in previous MC simulations of III- 
V DMS ESE1E1- We take L = 250-350a o , L w = 8a Q , 
and Ns — 0.8 — 1.6 x 10 5 Mn spins, which result in the 
desired Mn content x = 0.04. Hole concentrations in the 
range p = 0.1 — 1.4 x 10 11 cm~ 2 are experimentally rele- 
vant and correspond to Nh = 5—37 carriers in the simula- 
tion cell 0]. Periodic boundary conditions are employed 
in the in-plane directions and an infinite potential well 
is assumed normal to plane. For a given configuration 
of the spins, the Hamiltonian Hke is diagonalized in a 
plane-wave basis with two-dimensional (2D) wavevectors 
up to the truncation radius k c = 3.5(2ir/L) I 7|. The 
energy of the holes is determined by summing up the 
lowest Nh eigenvalues of H KE , thus neglecting carrier 
correlations and thermal excitations. The total energy is 
then obtained by adding the energy of the Mn-Mn AF 
interactions. Thermal fluctuations are simulated by the 
perturbative Monte Carlo method 0], with 2x 10 3 initial 
Monte Carlo sweeps (MCS) used for relaxation followed 
by 10 4 MCS for thermal averages. 
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FIG. 2: (a) Inverse Mn susceptibility of an undoped (Hke = 
0) Cdo.96Mno.04Te quantum well. Open symbols : Jij = and 
closed symbols : Jij 7^ 0. The intercept of the MC data and 
the horizontal (dotted) line at 1/xs = 1-83 (Eq.Q)gives the 
mean-field value of the critical temperature (T c ) for the doped 
system, (b) Variation of T c with doping level as obtained by 
MC simulations. Horizontal lines are the mean-field results, 
as obtained from (a). Dashed line : Jij = 0, solid line : 
Jij + 0. 



We show in Fig. 1 the temperature dependence of the 
spin polarization (|cr z |) = (| Ei°f l/-^0 and the suscepti- 
bility x* = [(I E* <l 2 > " (I <l) 2 ] /NT, where o = s 
or S and N — Nh or Ns, of a Cdo.96Mno.04Te QW at zero 
field and various hole concentrations. Non-zero values of 
spin polarization at high temperatures, of the order of 
l/2y/Nh, are due to finite size effects in carrier statis- 
tics. The abrupt increase in the hole and Mn spin polar- 
izations and the peak in the susceptibilities at the same 
temperature indicate a second order phase transition to a 
low-temperature state. While the hole spins become sat- 
urated immediately below T c , the much denser Mn spin 
subsystem saturates rather slowly on lowering tempera- 
ture. Notice that the range p = 0.10 — 1.40 x 10 11 cm" 2 
extends from the weakly to the strongly localization 
regime, as evidenced by the spatial extent of the carrier 
wave functions, discussed below. Interestingly, sponta- 
neous hole polarizations are observed in both these lo- 
calization regimes. 

Our MC data make it possible to verify the validity of 
the mean-field approximation (MFA) that was previously 
employed to describe the experimental findings 0,0]. In 
the absence of hole scattering (disorder), the MF T c is 
determined by the condition [2| 



A F N oX /3 2 p{e F )/4k B = l/xs(T). 



(1) 



Here Ap is the Fermi liquid parameter; p(s F ) is the den- 
sity of the carrier states at the Fermi level, where for 
the 2D case p{ep) — m* J dz\ip(z)\ 4 / nh 2 , and Xs{T) 
is the Mn spin susceptibility in the absence of carriers. 
Thus, within the MFA the ferromagnetic order appears 
if Xs^T) < 1.83 K, as obtained from Eq.JTJ for the 
sample in question and noninteracting holes (Ap = 1). 
Fig. 2 presents Xs(T) as determined by MC simulations 
for an undoped sample (Nh = 0), and implies that in 
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FIG. 3: Upper branch of the magnetization hysteresis loop of 
holes and Mn ions in a Cdo.96Mno.04Te quantum well at T = 
0.7T C . The hole density is p = 0.41 x 10 11 cm" 2 (iV h = 21). 

the relevant temperature range the AF interactions re- 
duce Xs(T) almost threefold from the value S 2 /3T of 
the paramagnetic case, in agreement with experimen- 
tal results for intrinsic Cdo.96Mno.04Te 0, I12I 11.^. As 
shown in Fig. 2(a), the Xs(T) data predict a MFA value 
of T c for the doped system equal to 1.25 ± 0.1 K and 
3.81 ±0.1 K when the AF interactions are switched on 
and off, respectively. A comparison of these values to 
the MC results (Fig. 2(b)) demonstrates that the MFA 
overestimates T c by a factor close to two if there are no 
AF interactions. This could be expected given the finite 
range of the RKKY-type interactions involved. However, 
in the presence of the AF interactions the difference be- 
tween the MC and MFA values of T c is much reduced. 
This implies that the system assumes the low tempera- 
ture phase at significantly higher temperature than that 
expected for a simple ferromagnetic order. 

To elucidate further the role of AF interactions we note 
first that the competition between long-range ferromag- 
netic interactions and short-range AF interactions was 
suggested to result in the formation of 180° domains, 
characterized by a wide spectrum of size distribution [Tlj | . 
This suggestion was put forward in order to explain the 
absence of magnetic hysteresis and, thus, of spontaneous 
global magnetization in p-(Cd,Mn)Te quantum wells be- 
low T c . In order to investigate this conjecture we have 
performed MC simulations of the field dependent mag- 
netization, shown in Fig. 3. When short-range AF in- 
teractions are switched off, the global spin polarization 
exhibits a square hysteresis specific to the Ising model 
in a magnetic field. However, in the presence of the AF 
interactions, a suppression of the remanent magnetiza- 
tion occurs and a relatively high value of the saturation 
field H sat ~ 150 Oe is obtained, in agreement with the 
experimental findings [lollllT|. Our results indicate that 
the AF interactions play a role of random fields that are 
known to turn the hysteretic behavior of Ising spins into 
a linear and reversible dependence of magnetization on 
the external field j^. 

An important aspect of our simulations is that they 
can provide information on the actual evolution of carrier 
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FIG. 4: Snapshot (at the final MCS) of the in-plane hole 
wave functions for the lowest eighteen eigenstates in a 
Cdo.96Mno.04Te quantum well with p = 0.41 x 10 11 cm -2 
(N h = 21) at high temperature (T = 5T C ). AF interactions 
are included. The grid scale is equal to 3.5a . 



wave functions with hole density, temperature and time. 
We plot in Fig. 4 the in-plane charge distribution of the 
carrier wave functions in a Cdo.96Mno.04Te quantum well 
with hole density p = 0.41 x 10 11 cm' 2 (N h = 21). These 
snapshots correspond to the final MC sweep of simula- 
tions carried out at high temperature (T = 5T C ) and 
zero field. The first few (n ~ 10) states exhibit strong 
in-plane localization, with a radius R ~ 20 — 30 nm, while 
higher states are weakly localized. Given that the one- 
hole states in our system are only weakly dependent on 
the total carrier density, we deduce that in lightly doped 
samples (p < 0.19 x 10 11 cm~ 2 , or Nh % 10) all states are 
strongly localized while in denser samples both strongly 
and weakly localized states exist. 

Furthermore, while the character of the wave functions 
varies appreciably with carrier density, it does not change 
substantially across T c . In particular, the sequence of 
strongly and weakly localized states, observed above T c 
(Fig. 4) remains unchanged below T c . This observation is 
consistent with the holes being scattered off the magnetic 
moments of Mn ions, whose ordering proceeds rather 
slowly on lowering temperature (Fig. 1), causing a corre- 
sponding slow variation of the scattering potential land- 
scape. The above features compose a carrier-mediated 
(RKKY-type) ferromagnetic transition in a system with 
disorder-driven Anderson localization. As regards the 
time evolution of the charge distribution, at T <C T c the 
hole gas remains fully spin polarized, while at T < T G our 
simulations indicated a slow and almost coherent rever- 
sal of the hole polarization. In this temperature regime 
the splitting of the PL line, being proportional to the 
absolute carrier polarization, will be non-zero, while its 
degree of circular polarization will average to zero. This 
behavior has indeed been observed in microluminescence 
studies (riT |. 

In order to assess possible effects of disorder we note 
that since in a true 2D system the density of states does 
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not vary with the Fermi energy, the MFA values of T G 
are expected to be independent of the hole concentra- 
tion. However, a lowering of T c when reducing p has 
been detected experimentally in thin QWs and attributed 
to disorder scattering |l2j. This behavior is reproduced 
by our simulations in the low-p limit (Fig. 2(b)), while 
the disorder-induced localization of the holes in the low- 
p limit (Fig. 4) is consistent with the conjecture on the 
role of disorder on the lowering of T c . At the same time, 
the lack of a substantial variation of the wave functions 
with temperature implies that the formation of magnetic 
polarons does not account for the phase transition in 
question. This conclusion supports the notion that the 
ferromagnetic instability driven by the RKKY-type in- 
teractions precedes the self-trapping of holes and the 
formation of magnetic polarons @- Finally, we note that 
the static spin susceptibility of the 2D carrier gas does not 
depend on the wave vector q, which suggests a tendency 
towards a spin-Peierls instability in such systems, occur- 
ring at q = 2fcp. Actually, such a spin-density wave sce- 
nario has been put forward in order to explain the above- 
mentioned absence of hysteresis in p-(Cd,Mn)Te quan- 
tum wells below T c If this would be the case, the 
spatial modulation of the carrier wave function should de- 
pend on temperature and, furthermore, its period should 
decrease as the square root of the carrier density. In 
contrast, our simulations point to an increase of the lo- 
calization radius with the hole concentration. 

In conclusion, we have simulated a quasi-2D system 
of carriers coupled to fluctuating localized spins in the 
regime corresponding to carrier-driven ferromagnetic in- 
stability and Anderson localization. Our results demon- 
strate that the RKKY-type interaction gives the domi- 
nant contribution to spin ordering, while the formation 
of magnetic polarons or spin density- waves appear to be 
of lesser importance. Furthermore, our findings iden- 
tify limits of the validity of the MFA in describing the 
hole-mediated ferromagnetism, particularly in the pres- 
ence of competing short-range AF interactions. Accord- 
ing to our results, the AF interactions decrease T G less 
than expected within the MFA but reduce strongly the 
remanence and the coercive field, which explains hitherto 
puzzling experimental results for p-(Cd,Mn)Te QWs. In- 
terestingly, these effects differentiate the (II,Mn)VI from 
the (III,Mn)V compounds, as in the former the ionized 
Mn acceptors attract holes and the AF interactions, even 
for first neighbors, become overcompensated by the fer- 
romagnetic coupling. Finally, it should be noted that 
carrier localization and spin-density wave instability in 
real systems are affected by the carrier-carrier repulsion, 
an effect neglected in the present simulations. However, 



we do not expect that the presence of carrier correlation 
will qualitative alter our conclusions. 
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